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ABSTRACT 

We re-analyse the ASCA and GINGA X~ray data from BY Cam, a slightly asyn- 
chronous magnetic accreting white dwarf. The spectra are strongly affected by complex 
absorption, which we model as a continuous (power law) distribution of covering frac- 
tion and column of neutral material. This absorption causes a smooth hardening of 
the spectrum below ^ 3 keV, and is probably produced by material in the preshock 
column which overlies the X~ray emission region. The ASCA data show that the in- 
trinsic emission from the shock is not consistent with a single temperature plasma. 
Significant iron L emission co-existing with iron K shell lines from H and He-like 
iron clearly shows that there is a wide range of temperatures present, as expected 
from a cooling shock structure. The GINGA data give the best constraints on the 



maximum temperature emission in the shocked plasma, with kT„ 



= 21 



+ 18 



keV. 



Cyclotron cooling should also be important, which supresses the highest temperature 
bremsstrahlung components, so the X-ray data only give a lower limit on the mass of 
the white dwarf of M > Q.5Mq. Reflection of the multi-temperature bremsstrahlung 
emission from the white dwarf surface is also significantly detected. 

We stress the importance of modelling all these effects in order to gain a physically 
self-consistent picture of the X-ray spectra from polars in general and BY Cam in 
particular. 
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1 INTRODUCTION 

Magnetic cataclysmic variables (polars or AM Her stars) are 
binary systems where a magnetised {B ~ lO'^ G) white dwarf 
accretes from a low mass companion (see e.g. the review by 
Cropper 1990). Such magnetic fields are strong enough to 
disrupt disk formation. Instead, the accreting stream is en- 
trained by the magnetic field and falls freely through the 
gravitational potential until it hits the white dwarf surface. 
The resultant strong shock has a typical temperature of tens 
of keV for optically thin material, giving rise to an X-ray 
emitting plasma. Normally such systems are locked into syn- 
chronous rotation by the magnetic field, so that (to zeroth 
order) the stream impacts onto the same part of the mag- 
netosphere. However, there are now three systems known 
where the orbital period and white dwarf spin period are 
slightly different, BY Cam, V1500 Cyg and RXJ 1940.2- 



1025 (e.g. Watson et al 1995), so that the accretion geometry 
changes continuously. 

Recent progress in understanding the detailed shape of 
the X-ray emission has concentrated on the effects of repro- 
cessing of the hard X-ray emission. Inclusion of the effect of 
Compton refiection from the white dwarf surface and com- 
plex absorption from the accreting material and the multi- 
temperature shock structure are necessary ingredients in the 
spectral model. While the importance of any one or two 
of these individual effects have long been recognised (e.g. 
Imamura & Durisen 1983 - absorption and shock structure; 
Swank, Fabian & Ross 1984; Beardmore et al 1995 - absorp- 
tion and reflection) all three are necessary in order to obtain 
a self-consistent description of the observed X-ray emission 
(Done et al 1995; Cropper et al 1997). Both the complex 
absorption and Compton refiection harden the spectrum. 
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leading to an overestimate of tlie plasma temperature. The 
line emission from the hot plasma gives another indepen- 
dent measure of the plasma temperature, from the ratio of 
6.7 and 6.9 keV lines from He- and H-like iron, respectively. 
The temperature of the line emitting material is then lower 
than that inferred from the 'observed' continuum. This dis- 
crepancy is exacerbated in detectors of low energy resolution 
as Compton reflection produces a fluorescence iron line at 
6.4 keV from the (nearly neutral) X-ray illuminated white 
dwarf surface, which is then blended by the instrument with 
the 6.7 and 6.9 keV lines. This reduces the mean iron line 
energy and leads to an underestimate of the plasma temper- 
ature. This systematic mismatch between 'observed' line and 
continuum temperatures is particularly evident in the large 
sample of GINGA (low spectral resolution) X-ray spectra 
of polars (Ishida 1991), even after the inclusion of complex 
absorption. Both complex absorption and Compton reflec- 
tion are required to obtain consistent temperatures from the 
line and continuum data (Beardmore et al 1995; Done et al 
1995), but then the derived temperature is of order 10 keV, 
much lower than the expected Rankine-Hugoniot tempera- 
ture for a strong shock of ~ 57{M/Mq){R/5.57x 10*)"^ keV 
(e.g. Frank, King & Raine 1992, Imamura et al 1997). Us- 
ing multi-temperature shock models rather than assuming a 
single temperature plasma goes some way to resolving this, 
as the mean temperature of a ~ 30 keV multi-temperature 
shock structure is ~ 10 keV (Imamura & Durisen 1983; Done 
et al 1995; Wu, Chanagen & Shaviv 1995; Cropper et al 
1997). Cyclotron cooling of the hottest /least dense material 
probably accounts for any remaining temperature discrep- 
ancy (Wu et al., 1995; Woelk & Beuermann 1996). 

While much real progress has been made, problems still 
remain. The most pressing of these have to do with the pre- 
shock accretion column. This material is irradiated by the 
X-ray emission, so should have a complex ionization struc- 
ture (see e.g. Ross & Fabian 1980; Swank et al 1984; Kallman 
et al 1993, 1996). It is also spatially extended over the X- 
ray source so that different segments of the X-ray emission 
travel through different path lengths of the material (Done 
et al 1995), especially as the accretion column is probably 
arc-like in cross-section, rather than circular (e.g. Cropper 
1990). The column should also add further complexity from 
secondary emission (Kallman et al 1993, 1996) and scatter- 
ing (Ishida et al 1997). Models including all these effects are 
currently being developed (Rainger et al 1997, in prepara- 
tion), but the situation is clearly complex. Even more de- 
grees of freedom are possible in the models: the material in 
the column is probably highly inhomogeneous, as dense 'rib- 
bons' of material embedded in a much less dense medium 
are required in order to explain the soft X-ray excess seen in 
several of these objects (e.g. Kuijpers & Pringle 1982; Frank, 
King & Lasota 1988; Ramsay et al., 1994). 

BY Cam has been extensively studied by previous 
X-ray instruments, and its spectrum shows many of the 
temperature discrepancies discussed above. Kallman et al 
(1993, 1996) propose an explanation for these in terms of 
secondary emission from the column, producing an addi- 
tional source of ionised 6.7 keV iron line emission. Here 
we reanalyse the GINGA and ASCA data from BY Cam 
and show that the spectrum is well described by the 
absorption/reflection/multi-temperature shock picture de- 



veloped above, and that the contribution of the pre-shock 
material to the iron K line emission is probably negligible. 



2 OBSERVATIONS 

The GINGA data were taken in 1988 February 7-10, with 
background observations taken directly before and after 
this. However, the second background could not be used as 
the ROSAT All Sky Survey Bright Source Catalog gives 3 
sources within 35 arcmin of the pointing position with total 
flux 1/5 of that from BY Cam itself. The first background 
alone is not sufficiently long to sample all the background 
conditions seen during the source observation so we exclude 
those parts of the source data which are not well matched by 
the background observation. Standard selection and clean- 
ing criteria were then applied to both source and background 
data, and the data were attitude corrected after a 'local' 
background subtraction (Hayashida et al 1989; Williams et 
al 1992), resulting in 46 ks of data. 

Ishida et al (1991) describes the time variability of the 
source during these observations. Clear (though generally 
partial) eclipses are seen in the first half of the dataset, but 
are absent in the second half. This is explained as a transi- 
tion in the accretion geometry from a two pole state, where 
only one of the accretion regions is periodically eclipsed by 
the white dwarf surface, to a state where the accretion orfly 
takes place onto the non-eclipsing pole (see also Silber et 
al 1992). We split the data into 3 spectra as in Ishida et al 
(1991), namely pulse high and pulse low (double pole ac- 
cretion with the eclipsing pole oriented towards and away 
from the observer respectively) and flaring (accretion onto 
a single non-eclipsing pole). 

The ASCA data taken on 1994 March 10-11 observa- 
tions and were previously reported by Kallman et al. (1996). 
FTOOLS 3.6 was used to generate standard response and 
arf files. The SIS event files were screened in bright2 mode 
with 1024-binning. The GIS3 data had to be binned onto 
128 channels as the data were taken during the time when 
the onboard CPU mode was wrongly set. This resulted in a 
total exposure time of 22000 seconds and total count rate of 
3.4 counts per second over the four detectors after the stan- 
dard reduction. The data were accumulated into one single 
spectrum as Kallman et al (1996) show no obvious hardness 
ratio changes. We restricted the energy range to 0.6-10 keV 
for SIS data and 0.8-10 keV for CIS data, and rebinned to 
have at least 20 counts per bin as required for validity of the 
X^ statistic. The relative normalization of the detectors was 
within 10 % from our preliminary analysis, thus we fitted 
all the spectra as a single data group. 

Kallman et al (1996) show the data folded on the po- 
larization period of Mason, Liebert & Schmidt (1989) with 
arbitrary phase. The orbital coverage is ~ 85% complete, 
and shows no evidence for an eclipse. The gaps in phase 
coverage are in two segments, each of width ~ 0.07 — 0.08 in 
phase. The GINGA eclipse in the pulse state has width > 0.1 
in phase (Ishida et al 1991), so it seems unlikely that the 
ASCA data could have avoided the eclipse altogether. The 
lack of any convincing ingress/egress in the ASCA data then 
support (although they cannot confirm) the view that the 
accretion geometry in BY Cam at the time of the ASCA ob- 
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servation was analogous to that of the GINGA flaring state 
data. 



3 SPECTRAL FITTING 

Version 9.1 of XSPEC (Arnaud 1996) is used to fit the spec- 
tra. The mekal code is used to model a single temperature 
hot plasma, while a modified version of the cevmkl code (see 
Done & Osborne 1997) is used for the multi-temperature 
plasma models which uses Morrison & McCammon (1983) 
abundance ratios. This gives a spectrum 



/M 



lO^-^K 



e{n,T) VT, 



T 



dT 



where e(n, T, u) is the bremstrahlung spectral emissivity at a 
given density n and temperature T while e(n, T) is the total 
(frequency integrated) emissivity (see also Done et al 1995). 
Models with a = 1 give a spectrum where each temperature 
component is weighted by its cooling time, i.e. the spectrum 
expected from pure bremsstrahlung cooling at constant pres- 
sure and gravitational field (see the discussion in Done et al 
1995 and references therein). 

The reflection code is described in Done et al (1995, 
1997), and here we assume a fixed inclination angle to the 
refiector of 60°. Error ranges are given as Ax^ = 2.7 unless 
otherwise stated, although this is an underestimate of the 
true 90% confidence limits where parameters are correlated. 

3.1 ASCA data 

Single temperature plasma models, where the abundances 
of all the elements heavier than He are assumed to scale 
together, give fairly good fit, with x^ = 1458/1272, but 
the minimum temperature of kT = 2OOI35 ^^^ shows that 
the continuum is much harder than can plausibly be ex- 
pected for optically thin accretion onto a white dwarf. In- 
clusion of an a = 1 multi-temperature continuum model 
does not significantly improve the fit (x^ = 1455/1272), 
and the derived (maximum) temperature is even higher. 



with kTmax 



940: 



keV, as the multi-temperature spec- 



trum includes the softer cooling components. Clearly there 
is some distortion present that is hardening the observed 
spectrum in the ASCA bandpass. Absorption is the most 
obvious way to do this since Compton reflection only con- 
tributes significantly to the spectrum above 5 keV, so is un- 
likely to strongly affect the ASCA data. Simple absorption 
(complete covering by neutral material) is already included 
in the fit, so we try partial covering by neutral material. 
Ishida (1991) shows that such complex absorption is not 
uncommon in polars, and uses this model to fit the GINGA 
data of BY Cam (Ishida et al., 1991). This gives a signif- 
icantly better fit for the multi-temperature a — 1 plasma 
model, with xl = 1414/1270 for Nh = 2.3lo;8 x 10^^ cm*^ 
and Cf = 0.20 ± 0.05 but still results in a physically un- 
reasonable maximum shock temperature of kT — 165lgg 
keV. A more realistic shock temperature of kT — 63l jg keV 
can be obtained with the single temperature plasma model 
with Nh = 2.7+1:0 X 10^^ cm~^ and Cf = 0.15 ± 0.03, but 
at the expense of physical consistency (we expect multi- 
temperature components to be present) and goodness of fit 
ixl = 1436/1270). 



The ASCA data extend to substantially lower energies 
than can be studied with GINGA, and so are more sensitive 
to the form of the absorber. A partial covering model is an 
ad hoc description of what is expected to be a much more 
complex situation. A power law distribution of column with 
covering fraction i.e. Cf{NH) oc Nf^ (see also Norton, Wat- 
son & King 1991) gives a good approximation to the complex 
absorption expected through neutral material extended over 
the source (e.g. Done et al 1995). The sum of covering frac- 
tions must be unity so that the emergent spectrum S{E) is 
related to the intrinsic spectrum Sint{E) by 

/•AfH.max 

SiE) = Sint{E)A N^exp{-NH(T{E))dNH (l) 

where 
A 



/3+1 



m 



73+1 

H,niax 



N, 



■13+1 
H,inin 



log(ArH,max/A'"H,min) 



...(3- 



and (t[E) is the photo-electric absorption cross-section from 
Morrison and McCammon (1983) abundances as given by 
the wabs model in XSPEC. We fix iYn.min = 10^° cm~^ 
and Figure 1 shows the transmitted fraction as a function 
of energy for Nn.max ~ 10'^'* cm~^ with /3 ranging from —1 
to 1. For comparison we also show the transmitted fraction 
obtained by the circular column model of Done et al (1995) 
for a mid column of 5 x 10^^ cm^^ (and hence a maximum 
column of 10^'* cm^^), and that obtained from complete cov- 
ering by a column of 5 x 10^^ cm~'^. The power law absorp- 
tion model gives a smooth hardening of the spectrum which 
matches the data significantly better {xt ~ 1399/1270) than 
the partial covering model. The derived maximum tempera- 
ture of a multi-temperature plasma is then the more physi- 
cally reasonable kTmax = 451^4 keV for a maximum column 
of 2.5l^g X lO^"* cm-^ with index /3 = -1.07±0.08. Figure 
2 shows the data and residuals for this model, along with 
the intrinsic spectrum before the absorption. 

Alternatively, the material may be ionized due to the 
intense irradiation of the pre-shock column by the X-ray 
emitting plasma. We use the XSPEC model absori to ap- 
proximate this, fixing the power law temperature of the illu- 
minated material to 3 x 10'' K and approximating the X-ray 
spectrum by a power law of photon index P = 1.5. Complete 
covering by a ionized material gives a worse description of 
the data, with xl = 1431/1270 for a column of 5.2+^;^ x 10^^ 
cm~'^ at an ionization parameter of f = 240l5^20i ^'^d again 
with an unphysically high temperature for the a — 1 plasma 
of kTmax ~ 270l|^QQ keV. Such models have been used above 
2 keV to match the complex absorption seen in the data 
(Bcardmore et al 1995; Cropper et al 1997), but they predict 
a large Oxygen edge which is inconsistent with the observed 
spectra below 1 keV. 

To summarise the results so far, the data strongly re- 
quire heavy absorption, which significantly hardens the spec- 
trum below 3 keV, but this absorption is complex. Complete 
covering of the source by a screen of material in a single ion- 
ization state (neutral or ionised) is strongly ruled out. In- 
stead the absorption can be described as a continuous range 
of columns through which the source is seen, as expected for 
the pre-shock material. While the ionization structure of the 
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Figure 1. The transmission as a function of energy through 
neutral material with a power law covering fraction - column 
distribution for /3 = —1,-0.5,0,0.5 and 1 from top to bottom. 
The dotted line shows the transmission through a circular neutral 
column where A^^ ^^.j = 5 X 10^^ cm~^ (Done et al 1995) while 
the dashed line shows the transmission through a complete screen 
with Nh = 5 X lO'^^ cm~^. The dash— dot line shows the transmis- 
sion through a circular column with radial density stratification 
(see Section 4.1). 



1 2 

channel energy (keV) 

Figure 2. The ASCA spectrum of BY Cam, modelled by a 
multi-temperature o = 1 plasma, absorbed by neutral material 
with a power law distribution of covering fraction with column. 
The top panel shows this model and the unfolded spectrum, to- 
gether with the intrinsic emission before absorption, while the 
bottom panel shows the remaining residuals to the fit. 



column is outside of the scope of this paper (see Rainger et 
al 1997), the data are clearly consistent with a power law 
column-covering fraction distribution of neutral material, so 
we use this model in all the following fits. The continuum is 
also significantly better described by a multi-temperature 
a = 1 cooling shock model than by a single temperature 
plasma, so we use this to describe the plasma emission in all 
the following fits. 

While the resultant fit is statistically adequate, there 
are clear residuals left around the iron K line (see Figure 
2). Adding a narrow Gaussian line gives a reduction in xt 
to 1363/1268 (significant at greater than 99.9 % confidence) 
for E = 6.39tJ',;'(5^ keV and equivalent width 86J:^^ eV. This 
additional line component is inconsistent with being at 6.7 
keV, as postulated by Kallman et al. (1996). A (nearly) neu- 
tral iron fluorescence line is expected from reflection of the 
intrinsic spectrum from the white dwarf surface. Fixing the 
line energy at 6.4 keV and including a neutral Compton re- 



flection continuum that must accompany any reflected line 
emission leads to Xv = 1360/1268, for an amount of reflec- 
tion R = l.Slj^'g (where R = 1 denotes the normalisation 
of the reflected continuuum expected from an isotropically 
illuminated slab covering a solid angle of 27r) assuming the 
abundances in the reflector are the same as those in the 
hot plasma, with maximum temperature kTmax = 38J;53 
keV. While the detection of the reflected continuum is only 
marginally signiflcant, its level is consistent with that ex- 
pected from the strength of the cold iron fluorescence line. 
The line equivalent width and amount of reflection contin- 
uum are anti-correlated in the fltting, so a self consistent 
reflection model (including both line and continuum) would 
be much better constrained. This flt is detailed in Table 1, 
and shown in Figure 3. 

The abundances of the elements in the hot plasma need 
not all scale together, especially as BY Cam is proposed to 
have abundance anomalies to explain the enhanced NV line 
emission observed in its UV spectrum (Bonnet-Bidaud & 
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Figure 3. The ASCA spectrum of BY Cam, modelled by a 
multi-temperature o = 1 plasma and its reflection from the white 
dwarf surface, absorbed by neutral material with a power law dis- 
tribution of covering fraction with column. The top panel shows 
this model and the unfolded spectrum, together with the intrin- 
sic emission before absorption, while the bottom panel shows the 
remaining residuals to the fit. 



Mouchet 1987). The N X-ray lines are below the 0.6 keV 
well calibrated low energy limit of the ASCA detectors, but 
can be investigated by including the less reliable data down 
to 0.4 keV. The soft excess is then also apparent (Kallman 
et al 1996), so a blackbody is needed to model this. The de- 
rived limit on the N abundance is then < 15x solar, which 
is not restrictive. Going back to the standard bandpass, the 
strongest lines expected at these plasma temperatures are 
Fe, Si, S and O. Letting the abundances of all these be free 



gives Af, = 0.41+^;«, As^ = 0.691^;^^, As = 0.39tl%, 
Ao = 0+°-^'^ and Arest = O.SSl^jg with xl = 1353/1264, 
i.e. there is no significant improvement in the fit as the in- 
clusion of 4 extra free parameters gives a reduction of only 
Ax^ = 7. Thus there is no strong evidence for non-solar 
abundance ratios of the elements, although O is marginally 
lower than expected. 

The intrinsic multi-temperature plasma emission can 
also be investigated in more detail. We have so far assumed 
that the plasma cools from the shock temperature to the 



white dwarf photosphere at ~ lO^'"" K ~ 27 eV, and that all 
these cooler components contribute to the spectrum. How- 
ever, this is not necessarily the case since the cooling plasma 
may be dense enough to become optically thick, so giving 
an apparent minimum temperature to the cooling radiation 
which is rather higher than that of the white dwarf. We 
modify the plasma emission model to include the minimum 
temperature as a free parameter of the fit, and find that this 
converges to an identical fit as before, with kTmin = 27 eV. 
The limits on the O abundance are also unchanged, with 
Ao = 0+°-^^ (xl = 1353/1266), showing that incomplete 
cooling is not responsible for any potential deficit of O line 
emission. 

Similarly, the temperature distribution need not be 
given by the a = 1 power law model expected from pure 
X-ray line and continuum cooling as the assumptions of 
bremsstrahlung only cooling at constant pressure and grav- 
itational field may not be accurate. Letting a be free gives 
a significant decrease in Xv to 1354/1267 for a = 0.6 and 

'^J-max -— <^<^ Rev. 

Thus the data are consistent with multi-temperature 
emission from complete cooling behind the shock, where the 
accreting material is ~ 0.5 x solar abundance in all the ele- 
ments. This intrinsic spectrum is then modified by refiection 
from the white dwarf surface, producing both reflected con- 
tinuum and associated 6.4 keV iron fluorescence line, and 
further distorted by complex absorption. The intrinsic emis- 
sion is required to be a multi-temperature plasma, and the 
cooling components are significantly detected in the data. 
Parameters for both the expected a = 1 cooling model and 
the better fitting model where a is free are given in Table 1. 
The best fit single temperature model (including refiection 
and complex absorption) is also tabulated for comparison, 
but is a worse fit by Ax^ = 20 as the data contain signifi- 
cant iron L line emission which cannot be fit by the single 
temperature models. This is the first observational confir- 
mation of the theoretically expected cooling of the shocked 
plasma in polars, although the distribution of cool compo- 
nents with temperature is marginally inconsistent with the 
predicted a = 1. 

3.2 GINGA data 

All 3 GINGA spectra give a lower fiux than the ASCA data, 
even in the overlapping 2-10 keV band. Ratios of these spec- 
tra with the best fit ASCA model are given in Figure 4. All 
the spectra are also systematically different in shape com- 
pared to the ASCA data, being generally softer. Thus there 
are significant spectral as well as intensity changes between 
the GINGA and ASCA observations. 

We first fit the GINGA low and high pulse state data, 
where the inclination angle to the eclipsing spot is changing 
with phase. We fit the two datasets simultaneously, con- 
straining the interstellar absorption, plasma temperature 
and abundances to be the same. Ishida et al (1991) show 
that complex absorption is required by the data, and model 
this by partial covering. A single temperature plasma model 
with partial covering gives x^ ~ 37/51. The low state spec- 
trum is significantly softer than the high state at low ener- 
gies, and in fact does not require any absorption. Constrain- 
ing the two partial coverers to the be same gives a much 
worse fit with x^ = 56/53. Replacing the partial coverer by 
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Table 1. ASCA spectral fitting with multi-temperature (pIT) and single temperature (IT) models 



Model 



N, 



H.Gal 



kT 



(keV) 



A' 



N 



xlO- 



R 



EW (eV) Af, 



H,niajc 



xt 



plT 
plT 

IT 



0.073 



7+0.03 
-0.04 

0.09 ±0.04 
0.06 ±0.03 



38 
55 
18 



+34 
13 

+ CX3 

-22 

+5 

6 



0.52 
0.45 
0.49 



+0.33 
0.15 

+0.30 
0.20 
0.20 
0.18 



7.0 
5.1 

3.5 



10 

0.60 ±0.25 



1 c+o.e 

J^-3-1.3 
-^■"-0.6 

1.1 ±0.6 



75 ±25 
75 ±25 
85 ±25 



5 

4.61 
1000 



2.6 

+7.0 

2.6 



a+0.12 

-0.13 

",+0.25 

-0.15 

-1.23 ±0.09 



-1.08: 

-1.0-t 



1360/1268 
1354/1267 
1380/1268 



" This is the maximum temperature in the multi-temperature models 

^ with respect to Morrison and McCammon (1983) 

'^ The normalisation of the plasma model in units of 10~^''/(47rD^)x Emission Measure 

<* Units of 10^2 (.m-2 

" Parameter fixed 



ma. = 5.0+^:^ and 770+^:^0 x 10^ 
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Figure 4. The ratio of the GINGA spectra to the best fit q = 1 
ASCA spectral model shown in Figure 3. The low pulse phase, 
high pulse phase and flare spectra are shown in the top, middle 
and bottom panels, respectively. 



the power law column-covering fraction model used for the 
ASCA data gives an equally good fit, showing that the wider 
ASCA bandpass is needed to properly constrain the complex 
absorption properties. We use the power law absorber for 
consistency with the ASCA results, but since A'^n.max and /3 
are strongly correlated and poorly determined we constrain 
(3 to be equal across all the datasets. This gives x^ ~ 35/52 



for /3 = -1.01 ± 0.08, iVH,ma:r = 5.0l^:g and 770 
cm~^ for the low and high state spectra, respectively. 

Again the residuals clearly indicate structure around 
the iron line and edge energies so we fit a reflection contin- 
uum spectrum and 6.4 keV iron fluorescence line. This gives 
a significant improvement in the fit with x? ~ 28/50. The 
continuum was then replaced with the multi-temperature 
a = 1 emission model. This did not give a significantly 
better fit to the data, showing that the 2-20 keV GINGA 
bandpass is not sufficient to distinguish between single and 
multi-temperature emission models. However, we chose this 
model for further study, for ease of comparison with the 
ASCA fits. The fiaring state spectrum is then also fit with 
this model, and the results detailed in Table 2. Error ranges 
on the absorber are not given since this quantity is uncon- 
strained by the data. The harder spectrum seen in the high 
state can equally well be described by a larger contribution 
from the refiection continuum as by more absorption. 



4 DISCUSSION 

4.1 Absorption 

The complex absorption is clearly a major component in 
the spectrum, giving a smooth hardening at energies lower 
than ~ 4 keV which cannot be described by complete cov- 
ering by material of a single ionization parameter. The lack 
of the standard absorption signatures (strong low energy 
cutoff or ionised edges) means that it is not always obvi- 
ous that this component is present in the spectrum. In BY 
Cam its main signature is that the derived temperature for 
the ASCA spectrum becomes unphysically high without it. 
We caution that where this complex absorption is present, 
the observed bremsstrahlung luminosity can be severely un- 
derestimated in spectra which only extend to 2-3 keV (e.g. 
ROSAT; Ramsay et al., 1994). 

Our model for this absorption is one in which neu- 
tral material with a range of covering fractions obscures the 
source, so that Cf{N-ii) oc A^^. The simplest geometry of a 
cylindrical neutral column overlaying a circular source pre- 
dicts a distribution of column with covering fraction such 
that (3^0 (see Figure 1). This is very different from the 
observed /3 = — 1, indicating that the situation is (unsur- 
prisingly) much more complex. 

The absorption can only be constrained by the ASCA 
data which is phase averaged, so some of the discrepancy 
in expected column distribution could arise from co-adding 
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Table 2. Simultaneous spectral fits to the GINGA pulse low (Low), pulse high (High) and flaring data with a multi-temperature plasma 
continuum, reflection and complex absorption. The temperature of the High and Low pulse states are tied together, while the abundance, 
temperature distribution and absorption power law are tied across all the datasets. 



Data 



N 



H.Gal 



(keV) 



a[ 



N 



xlO" 



R EW (cV) AfH,max0 



/3 



Total xl 



'^■^^-0.09 



Low 

High 

Flaring 



0.0+0' 



2ltf 
25tf 



2.4 
5.2 
6.5 



1 1+1-2 
K+1.2 
-1.3 
1 0+1° 



1.6-* 



0+135 

35 



45 



+70 

35 

+80 



1.9^ 



-0.9 41.4/73 



4.0 



° This is the maximum temperature in the multi-temperature models 

^ with respect to Morrison and McCammon (1983) 

'^ The normalisation of the plasma model in units of 10^i*/(47rD2)x Emission Measure 

<* Units of 10^2 cm-2 

■^ Parameter fixed 



data. Further distortions could also be produced if the accre- 
tion column is ionised. This would have the effect of making 
the absorption from a circular column less dramatic than 
predicted from neutral material, as required by the (5 = —1 
distribution. Photo-ionization of a uniform accretion col- 
umn has long been known to be likely (Ross & Fabian 1980; 
Swank, Fabian & Ross 1984) but is much less probable if the 
column has a highly inhomogeneous 'blobby' structure (e.g. 
Frank, King & Lasota 1988). In this case radial structure 
can alter the observed properties. The observed P = —1 
distribution can be reproduced by a model in which neu- 
tral blobs accrete preferentially towards the centre of the 
column, with blob density oc r~^, so that the much less of 
the X-ray source is covered by the highest column density 
material (dash-dotted line in Figure 1). 

All accretion column models predict a substantial 
change in absorption with phase. For a neutral, constant 
density column, there should be much less absorption when 
the hard X-ray emission region is on the limb of the white 
dwarf (pulse low state) than when it is viewed more directly 
(pulse high) as the mean path length through the accretion 
column is smaller (Done et al 1995). Allowing for ionization 
enhances this effect, as the high ionization material (and 
hence weakest absorption) is produced in the layers closest 
to the shock. The GINGA data are indeed consistent with 
a change in absorption between low and high state, in the 
sense that the low state is less absorbed (see Table 2), but 
a complete model of the density and ionization structure of 
the photo-ionised column to determine both its absorption 
and emission properties is clearly needed before such ideas 
can be properly tested. 

The absorption is also presumably linked to the mass 
accretion rate, as there is more material in the column at 
high mass accretion rates. Thus for a neutral column we ex- 
pect that the absorption should correlate with X-ray lumi- 
nosity. By contrast, an ionised column may show absorption 
anti-correlated with luminosity as the bremsstrahlung in- 
tensity is ex n^ while the column density is only oc n. Again, 
proper modelling of the column is required in order to deter- 
mine the observational consequences of this, and better data 
are required to test whether the absorber is predominantly 
ionised or neutral. 



4.2 Continuum 

A single temperature plasma continuum is clearly ruled out 
by the ASCA data. Multi-temperature plasma is expected 
physically, as the hot shocked gas must cool as it emits X- 
ray radiation. The observational signature of this is that 
the spectra show significant emission both from the H and 
He-like Kq lines at 6.9 and 6.7 keV and iron L shell lines 
around 1 keV from much lower ionization species. For all 
of these different Fe ion states to co-exist requires a multi- 
temperature plasma. The observed maximum X-ray tem- 
perature of 21^4 keV (from the better constrained GINGA 
fits) is rather lower than the ~ 57 keV expected from a 
pure bremsstrahlung shock above a solar mass white dwarf. 
Taken at face value this would limit the mass of the white 
dwarf to 0.6to'i- However, BY Cam has a substantial mag- 
netic field of ~ 28 MG (Schwope 1996 revised from the 
Cropper et al 1989 value of 40.8 MG), so cyclotron cool- 
ing may be important. The high temperature/low density 
material close to the shock preferentially cools via cyclotron 



(emissivity oc T 



Wu et al 1995) while at lower tem- 



peratures/higher densities the bremsstrahlung (emissivity 
oc T^''^n^) becomes more important. This has the effect of 
reducing the maximum X-ray bremsstrahlung temperature 
below that expected from a simple strong shock (Wu et al 
1995; Woelk & Beuermann 1996). 

Composite cyclotron-bremsstrahlung cooling shocks 
also have another feature not seen in pure bremsstrahlung 
shocks, namely that the maximum X-ray temperature 
changes as a function of mass accretion rate. Higher mass 
accretion rates increase the density of the shocked ma- 
terial, and so shift the balance between cyclotron and 
bremsstrahlung cooling in favour of bremsstrahlung (Woelk 
& Beuermann 1996, Wu et al 1995). Thus we expect the 
maximum shock temperature as seen from the X-ray emis- 
sion to increase with increasing mass accretion rate. We use 
the code of Wu et al (1995) to calculate the temperature at 
which bremsstrahlung and cyclotron cooling are equal for 
_B = 30 MG, M — Mq and mass accretion rate of 3.14 and 
6.28 X lO^^^ g s~i (corresponding to luminosity of 0.75 and 
1.5 X 10^'' ergs s~i, respectively). This gives Tmax ~ 25 and 
35 keV. These are intruigingly close to the best fit derived 
maximum temperatures for the a = 1 GINGA and ASCA 
data of 20 and 38 keV, respectively, for a factor of ~ 2 differ- 
ence in mean 2-10 keV count rate. While the large error bars 
mean that the effect is not significant here, these tempera- 
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ture changes should be seen in high quality, broad bandpass 
spectra from SAX and/or XTE 

Cyclotron cooling has a strong and variable effect in 
suppressing the high temperature X-ray emission. These BY 
Cam data stress again the point made by Wu et al (1995) 
and Woelk & Beuermann (1996) that equating the max- 
imum bremsstrahlung temperature to the shock tempera- 
ture is not a good guide to the white dwarf mass. However, 
downstream from the shock, once the density and tempera- 
ture are such that bremsstrahlung dominates the cooling 
then the shock structure is again given by the standard 
bremsstrahlung a — 1 cooling. This can be seen from the 
results of the code of Wu et al (1995) and also theoretically, 
since when bremsstralung cooling dominates all the equa- 
tions are the same as for a bremsstralung only shock. The 
derived temperature distribution seen in the ASCA data is 
(marginally) significantly different from this, with a ~ 0.6 
i.e. more cool components than expected. This may be an ar- 
tifact of the ad hoc absorption model used (see above). If the 
real absorber includes complex contributions from ionised 
material then some structure should be present at the OVIII 
K and Fe L line and edge energies which may distort the de- 
rived continuum. Alternatively, this may be a real effect, 
showing that the shock structure deviates significantly from 
that expected, perhaps due to conduction and/or Compton 
cooling (Imamura & Durisen 1983; Imamura et al 1987). 

4.3 Abundances 

The X-ray line emission from BY Cam is weaker than ex- 
pected from solar abundance (Morrison and McCammon 
1983) coronal models. This has been seen systematically 
from many classes of objects, including stellar coronae (see 
e.g. Guedel et al 1997; Mewe et al 1997) as well as other Cat- 
aclysmic Variables (Done et al 1995; Done & Osborne 1997). 
While deviations from coronal equilibrium may be expected, 
none of these seem able to produce the observed weak lines 
(Done et al 1995, Done and Osborne 1997), so it seems more 
likely that this represents a true under-abundance with re- 
spect to solar. This seems more plausible now that optical 
studies suggest stars in the solar neighbourhood have con- 
siderable dispersion in metallicity, and that the mean [Fe/H] 
is 0.25-0.3 dex lower than that of the Sun (Edvardsson et al 
1993), but is still an uncomfortable conclusion. 

Strong UV lines from NV have led to a claim that the N 
abundance is enhanced (Bonnet-Bidaud &: Mouchet 1987). 
X-ray line strengths from H and He-like ions give much less 
ambiguous determinations of abundances than UV or optical 
photo-ionised lines, but the X-ray N lines are below the well 
calibrated energy limit of the detector and the derived limit 
of < 15 X solar is not restrictive. Whether or not N is indeed 
enhanced, the suggested explanation of a recent novae explo- 
sion as the cause (Bonnet-Bidaud & Mouchet 1987) seems 
unlikely since the lines arise from the accretion stream i.e. 
material from the secondary star. Even if the secondary did 
manage to capture some fraction of the nuclear processed 
material, the elemental abundance anomalies would become 
diluted by mixing with the un-reprocessed material in the 
stellar atmosphere. However, recent HST spectra of the sur- 
face layers of the accreting white dwarf in VW Hydri also 
seem to show abundance anomalies, with Al~ 15 x solar, 
and P~ 900 X solar (Sion et al 1997)! Again, these have 



been interpreted as being indicative of a recent novae, so we 
examine other possible X-ray signatures of nuclear repro- 
cessing in BY Cam. 

The core composition of the accreting white dwarf can 
either be CO or ONeMg. For a CO white dwarf, calculations 
by Kovetz & Prialnik (1997) show that while N is almost al- 
ways enhanced in a nova explosion, the O abundance can 
be almost unaffected or strongly depleted, while the heav- 
ier elements are not substantially enriched. This is consis- 
tent with the ASCA observations. Conversely, for an ONeMg 
white dwarf core, N/0< 1 requires that the abundances of 
P and S should also be strongly enhanced (Politano et al 
1995). This type of nova explosion can then be clearly ruled 
out since S is not overabundant (see Section 3.1). Thus if a 
recent novae has enriched the N abundance, it would have 
to be from accretion onto a CO core, although a more likely 
explanation for anomalous C/N ratios is that the secondary 
is evolved (Mouchet et al 1997). 

4.4 Comparison with previous work 

Physically motivated (as opposed to phenomolgical) descrip- 
tions of the spectrum of BY Cam are given by Kallman 
et al (1996), using the same ASCA data as studied here. 
However, their conclusions are rather different in that they 
propose a strong line emission component at 6.7 keV from 
the pre-shock column, rather than absorption. They re- 
quire this because their assumed continuum form (a 30 keV 
bremsstrahlung) produces copious 6.9 keV line from H-like 
iron, but has a low He-like iron ion fraction, so cannot pro- 
duce much of the observed 6.7 keV line. We show that this 
conclusion is obviated by using a continuum given by the 
(physically expected) multi-temperature plasma emission 
from a cooling shock, as this can produce both the 6.7 and 
6.9 keV lines. Kallman et al (1996) discount this possibility 
due to the weakness of the observed 6.9 keV line compared 
to that expected from lower temperature plasma with solar 
abundances. This objection is removed by allowing iron (and 
the other element abundances) to be sub-solar. Kallman et 
al (1996) also comment that with lower temperature plasma 
the observed continuum is too flat. In our analysis this is 
resolved by the inclusion of complex absorption from the 
column. We note that Kallman et al (1996) require a pre- 
shock column of order --^ 10^'' cm~^ overlaying the X-ray 
emission region. Thus their model is not self-consistent as 
it does not take into account the strong absorption that this 
predicts. Our model is also inconsistent in not including the 
emission from the X-ray illuminated column, but the sub- 
solar iron abundances means that any predicted iron K line 
emission is reduced by a factor 3 below that of Kallman et 
al (1996) i.e. < 25 eV equivalent width, which is smaller 
than the error bars. The photo-ionised column is instead 
expected to produce the bulk of its emission from contin- 
uum/recombination continuum/lines at energies below ~ 1 
keV, as noted by Kallman et al (1993) where a large neu- 
tral absorption column is required to suppress the strong 
predicted low energy emission. 

Kallman et al (1996) also speculate on the presence of a 
reflection component to produce the observed 6.4 keV line. 
Their derived limit of Q/2tt < 0.6 is again assuming solar 
abundances, without including the inclination effects, and 
using the results of early reflection calculations, which gave 
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EW ~ 210il/27r eV for a 20 keV bremsstrahlung illumina- 
tion. More recent Monte-Carlo results (e.g. George & Fabian 
1991, Matt, Perola & Piro 1991, Van Teesling, Kaastra & 
Heise 1996), scaled to a 20 keV bremsstrahlung continuum 
(e.g. Beardmore et al 1995) indicate that even for a face 
on, solar abundance slab the line emission is no more than 
~ 150r2/27r eV, while for a mean (phase averaged) viewing 
angle of 60° this reduces to ~ 110r2/27r eV, while abun- 
dances of ~ 0.4 X solar reduces it further to ~ 90f2/27r eV 
(George & Fabian 1991). Thus the 6.4 keV line strength is 
easily consistent with reflection from a shock just above the 
white dwarf surface i.e. with Q./2-k ~ 1. 



5 CONCLUSIONS 

The BY Cam data from ASCA and GINGA give a physically 
self consistent picture of an accretion column shock, cool- 
ing by both cyclotron and bremsstrahlung emission. This 
multi-temperature X-ray continuum illuminates the white 
dwarf surface, producing a reflection continuum and iron 
fluorescence line. Absorption from the pre-shock column 
strongly modifies the observed spectral form, and is a signif- 
icant source of uncertainty since it depends on whether the 
pre-shock column is uniform, or blobby (and hence ionised 
or nearly neutral), circular or arc-like in cross-section and 
whether there is radial density structure. We model this ab- 
sorption by neutral material where the covering fraction is 
a power law function of the column. This could represent a 
physical situation where the column is blobby (unionised), 
with more blobs accreting towards the centre of a circular 
column. However, it is more likely that this form merely 
gives a suitable approximation to a more complex situation, 
and we urge further theoretical modelling of the pre-shock 
flow. 

The multi-temperature emission and reflection model, 
together with the power law neutral absorption model will 
be made publically available in the next release of XSPEC. 
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